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ABSTRACT 

We perform a set of non-radiative cosmological simulations of a preheated intracluster medium in 
which the entropy of the gas was uniformly boosted at high redshift. The results of these simulations 
are used first to test the current analytic techniques of preheating via entropy input in the smooth 
accretion limit. When the unmodified profile is taken directly from simulations, we find that this model 
is in excellent agreement with the results of our simulations. This suggests that preheated efficiently 
smoothes the accreted gas, and therefore a shift in the unmodified profile is a good approximation 
even with a realistic accretion history. When we examine the simulation results in detail, we do not 
find strong evidence for entropy amplification, at least for the high-redshift preheating model adopted 
here. In the second section of the paper, we compare the results of the preheating simulations to 
recent observations. We show - in agreement with previous work - that for a reasonable amount of 
preheating, a satisfactory match can be found to the mass-temperature and luminosity-temperature 
relations. However - as noted by previous authors - we find that the entropy profiles of the simulated 
groups are much too flat compared to observations. In particular, while rich clusters converge on the 
adiabatic self-similar scaling at large radius, no single value of the entropy input during preheating 
can simultaneously reproduce both the core and outer entropy levels. As a result, we confirm that 
the simple preheating scenario for galaxy cluster formation, in which entropy is injected universally 
at high redshift, is inconsistent with observations. 



Subject headings: cosmology: theory - intergalactic medium 
clusters 



methods: numerical - X-rays: galaxies: 
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1. INTRODUCTION 

Clusters of galaxies are potentially powerful tools 
for precision cosmology. In ACDM and other mod- 
els of hierarchical structure formation, the mass func- 
tion of clusters and its evolut ion are very sensi- 
tive to cosmological parameter s (iBahcall fc CenI 
Haiman. M ohr. fc Holderl 120011: IBahcall fc Bodd ' 



1991 



2001 



Younger. B ahcall. fc Bodd 120051) . Recently, many au- 
thors have used the observed cluster abundance in 
a variety of wavelength ranges to constrain both the 
present day mean matter density and the am- 
plitude of_Jineaj_^uchiations on 8 Mpc scales 
e.g., Eeni^ fc ArnaudI [T99lt Ivima fc Liddld UMl 



Borgani et al.^'^2001a"; Bahcall ct al."2003': Picrpaol i et al 
2003: Shueckcr ct al . 2003; .Henry ,2004; Gla ddcrs et al 
20061 ). In the future, the first large-scale X-ray 



surveys will yield cat a logs of many thousands of 
clusters (jHaiman et al.l |2005() that will place very 
tight constraints on cosmological parameters, includ- 
ing the dark energy equat i on of state and it s evo- 
lution (iMaiumdar fc Mohd 120041 : IWang et all 120041 : 
lYounger et al.ll2006f) . 

With such a potentially promising dataset on the hori- 
zon, particular attention has been paid recently to ad- 
dressing potential systematic un certaint i es tha t may bias 
measurements of cosmology. As iHenrvl (|2004[ ) observed, 
the systematic uncertainty in determining erg using cur- 
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rent cluster measurements is driven by the calibration 
of the mass-temperature relation. If cosmological results 
from future cluster surveys are to be robust, accurate 
theoretical models of cluster formation and evolution, in- 
cluding identification and proper treatment of the domi- 
nant heating and cooling mechanisms of the intracluster 
gas, are required. 

The observational properties of the intracluster 
medium (ICM) are driven primarily by simple grav- 
itational collapse. If this were the only important 
physica, ! mech anism, clusters would scale self-similarly 
(jKaiseil [l986l ). However, it has been known for some 
time that the observed bulk properties of galaxy clus- 
ters - temperature, mass, lum inosity - do not con- 
form to the self - simila r model (lEvrard fc Henrvl Il991t 



Edge fc Stewart 1991"; 'David et al. 1993; Markevich 
1998 : Arnaud fc Evrard .1999; Reipric h fc Bo hringer 
2002f ). This effect is also seen in the luminosity- 



temperature scaling of groups (jHelsdon fc Ponmanll2000( l 
and the faintn ess of the unresolved ^ 1 keV X- 
ray background (|Penl fl99l IWu. Fabian, fc NulsenI [20011 : 
iBrvan fc VoitI l2001h . These departures from self- 
similarity are affected by non-gravitational baryonic 
physics; a complex interplay between various thermody- 
namic proc esses involving star formation and galaxy evo- 
lution (see iRosati. Borgani. fc Nori^l2002l : TVoitll2005L 
and references therein). 

Motivated in part by observati ons suggesting a uni- 
versal entropy floor for clusters ( Ponman et al.l Il999t 



iLlovd-Davies. Ponman. fc CannonI |2000| ). manv studies 
have proposed "preheating" of the gas in or der to 
expla in this departure f rom s e ff-similarity (e.g., [Kaiseij 
119911: lEvrard fc Henrvl 119911 : iNavarro. Frenk fc Whiti 



1991 



ICavaliere. Menci. fc Tozzl 



J1997I : 



Balogh. Babul, fc Pattodll999l : iPonman et al.lll999l ). In 
this scenario, the effects of energy input into the ICM 
from non-gravitational processes such as supernovae, 
star formation, and galactic winds are approximated by 
a high-redshift entropy modification that parti cularly 
affects low-mass systems ( Boweijll997l: pTo zzi & N ormanI 
[2OOII: iVoit fc Brvanl [200l iVoit et alJ 1 20021 i200a. Sub 



— , : IT -41^ 

sequent numerical (IBorgani et al.l I2001bt IKav et al 

2004t IBorgani et~ar 2005; Muanwong, Kav. fc Thomas 



2006) and analytic (Babul ct al. 2002_; Voit et al.ll2002l : 



McCarthy et al.l[200 3a. 2004) modeling has been success- 



ful in reproducing many of the observed cluster scaling 
relations, both in the X-ray and Sunyaev -Zel'dovich 
(SZ) effect (ISunyaev fc Zeldovich' 1972'. 198^. 

In particular. iVoit et a l. (2002, 2003) propose that a 
simple shift of the entropy profile may be a good approx- 
imation to the effects of preheating.^ However, since 
this results was obtained in the smooth accretion limit, 
and real accretion of preheated gas is likely to be some- 
what lumpy, it may not hold generally. This is because it 
is not clear a-priori how efficiently universal preheating 
will smooth the gas distribution prior to cluster collapse. 
In this work we compare the predictions of this shift 
model for preheating to the results of high-resolution 
cosmological simulations. Such a comparison tests this 
analytic prescription, given a realistic accretion history. 
And, it will motivate other model assumptions, such as 
the choice of boundary conditions, that are more difficult 
to derive from first principles. 

At the same time, recent observations of the entropy 
profiles of nearby clusters have suggested that preheating 
scenarios which invoke a uni versal entropy floor are incor- 
rect in detail jP onman. Sa nderson, fc Finoeuenovl 120031: 



Pratt fc Arnaudi i2005t iPratt. Arnaud. fc Pointecouteaul 
2006D . The same simulations we use to test analytic pre- 



scriptions offer the added benefit of a fully informed com- 
parison of the predictions of the preheating scenario to 
observations. 

This work is divided into three sections. In § [2l we 
summarize the simulation and cluster identification pro- 
cedure, in § [3l we compare the corresponding model 
predictions to the results of our simulations. Finally, in 
§ m we compare the results of our simulations to recent 
observations of nearby clusters. Throughout this work, 
"preheating" should be taken to refer specifically to the 
shift model of IVoit et al.l (|2(iol [2003h . 

2. THE SIMULATIONS 

2.1. Overview 

The simulations were performed with the 
Adaptive-Mesh-Refinernent code ENZO (BrvanI 119991 : 
iNorman fc B7^fl999l : lO'Shea et al]l20?)l assuming a 
flat ACDM cosmology. The cosmolo gical model is consis- 
tent with the WMAP Year 3 results (|Spergel et al.ll2006f) 
with (a^,f7A,f^b,/i,cr8) = (0.25,0.75,0.046,0.7,0.75). 
We ran two complementary sets of simulations. The 
first simulation was performed in a comoving periodic 

^ Throughout this paper we use the term preheating to refer 
to a constant entropy increase, which corresponds to a spatially 
constant e nergy input at early times, when the gas is smooth. 
IVoit et al.l ||2002) also discuss other forms of modifying the en- 
tropy profile, largely motivated by the effect of cooling and star 
formation, which we do not discuss in this paper. 



box with L = 100 h^^ Mpc on a side with N = 256^ 
dark matter particles and an equal number of grids, 
allowing for six levels of additional refinement that 
yielded a minimum cell size of ~ 6 kpc. This level 
of resolution was found (via a smaller set of higher 
resolution simulations) to be sufficient for reproducing 
the core entropy level of the smallest clusters we will 
investigate in this paper. The initial conditions were 
generated with the fitting fo rm of the linear d a rk ma tter 
power spectrum given by lEisenstein fc Hul (1999) at 
an initial redshift of z = 60 (for d etails regarding the 
initialization of such simulations see lBertschingedll998l . 
and references therein). 

Due to its relatively small size, this 100 /i^^ Mpc sim- 
ulation produced only a limited number of high mass 
clusters. In order to fill out the higher-mass end of the 
scaling relations, we performed a second set of simula- 
tions using the same cosmology and method of generat- 
ing initial conditions. First, a dark-matter only run was 
performed in a comoving periodic box L = 400 Mpc 
on a side with N = 128'^ dark matter particles. Then, 
the largest dark -matter halos were ide ntified using the 
HOP algorithm (jEisenstein fc IIuflll998f ). and the simu- 
lation was re-run with two nested grids: the first in the 
original L = 400 h'^ Mpc box with N = 128^ dark- 
matter particles and an equal number of grids, and the 
second in a high-resolution region, centered on one of the 
haloes, with L ~ 50 Mpc and the same number of 
particles and grids, using the same refinement technique 
as the larger simulations. This gave our high-resolution 
region equivalent mass and spatial resolution to the orig- 
inal runs. We repeated this process for each of the four 
largest dark-matter haloes. 

2.2. Preheating the ICM 

IVoit et all (|2003l ) have suggested that simple preheat- 
ing in the smooth accretion limit is well-approximated 
by a universal shift of the unmodified cluster entropy 
profile. This corresponds, in the context of a full cosmo- 
logical simulation, to raising the entropy level of all the 
gas in the universe at high-redshift, well in advance of 
cluster collapse, but within the time frame during which 
we expect to find galactic and stellar entropy input to 
the ICM. Therefore, a fixed amount of entropy i^o was 
added at 2; = 10 by increasing the thermal energy of each 
grid-point by 



fcAT = 



Pa 



He 



7-1 



(1) 



where pg is the baryon density at that grid-point. We 
assume an ideal gas consisting of a fully ionized H-He 
plasma with the cosmic helium mass fraction Y[jo = 0.25 
and 7 = 5/3. This was done for four cases: no pre- 
heating [Kq = 0) and — 78, 155, and 311 



-1/3 
70 



keV cm^ . The effects of radiative cooling and other 
additional sources of non-gravitational heating were ne- 
glected. Snapshots of the z = emission-weighted 2- 
D projection of the gas temperature in the same clus- 
ter from two of the simulations are shown in Figure [TJ 
This figure demonstrates visually that adding entropy 
smoothes the small-scale distribution, but leaves the 
larger scale structures largely intact. 



Fig. 1. — Emission weighted 2-D projection of tlie 2 = gas temperature for a cluster witii M200 ~ 7 X 10"'^'' h^^ Mq two of the 
simulations: Kq = (left) and 311 h.^^^^ keV cm^ (right). Note how the preheating entropy input is effectively a smoothing operation. 



2.3. Cluster Identification and Analysis 

At each redshift ou tput, clusters were identi fied using 
the HOP algorithm of lEisenstein fc HutI ()1998D . yielding 
clusters at z = up to A'hoQ ~ 3x lO" h:^^ Mq in the 100 
Mpc simulations, and clusters up to M200 8 x 10^'* 
li^Q Mq in the large box simulations, where Af2oo is 
the total mass enclosed within a spherical overdensity 
of 200 relative t o the critical density. Th ese were then 
analyzed using a lRavmond fc SmithI (|1977f ) cooling func- 
tion assuming a fully ionized H-He plasma with abun- 
dance Z = Q.ZZq. The cluster center for large clusters 
(-^200 ^ 4 X 10^^ Mq) was calculated using an iter- 
ative process of re-centering on the center of mass using 
successively smaller radii. For all others, it was taken 
as the densest dark-matter particle from the hop out- 
put. We found that the choice of centering technique 
had no more than a 15% effect on the core properties of 
the cluster, even at the highest masses. 

For each cluster, we determined the mass {Tgw) and 
emission (Te^,) weighted gas temperatures, and the bolo- 
metric X-ray luminosity {L{,oi)- In addition, we consider 
observables associated with the Sunyaev-Zel'dovich ef- 
fect (SZ; Sunyacv fc Zcldovich 1972, 1980), including the 
central decrement (j/c) and frequency independent, angu- 
lar diameter distance scaled (FIADS) integrated SZ lu- 
minosity out to i?2oo {Syd\/ fy). These were calculated 
according to 

= I Pede, (2) 



where err is the Thompson cross section, TOg is the mass 
of the electron, c, is the speed of light, = n^kT is the 
electron pressure, and di is the line element along the 
line of sight through the cluster center, and 



Si/d\ 



[kTcMB? 



(JT 



[he) 



PedV, 



(3) 



where k and h arc the Boltz mann and Planck co nstants 
respectively, Tcmb = 2.728 (Fixse n et all I1996D is the 
Cosmic Microwave Background (CMB) temperature, and 
the integral is performed over a spherical volume out to 



-^200- Finally, we calculated the gas-weighted entropy 
profile, where entropy was defined as 



K. 



gw 



'^^ gw "'e 



where Ue 



(4) 

'TiT' "^ 2^° average number density of 

electrons in a given shell, again assuming a fully ionized 
plasma. 

3. COMPARISON TO ANALYTIC MODELS 

3.1. The Analytic Model 

Our anal ytic model follows the techniques of lVoit et al.l 
(J200^ and lYounger et all ()2006f ). Given an initial, un- 
modified, spherically symmetric model in hydrostatic 
equilibrium, then entropy profile is "modified" - in this 
case with a simple shift at all radii - and the equations of 
hydrostatic equilibrium are re-integrated, yielding a final 
"modified" gas and temperature distribution. Therefore, 
if K represents the shifted entropy profile, both the un- 
modified and modified models satisfy 



dP 

— ^rjg{r)pg{P,K) 



(5) 



where, casting the pressure and temperature in terms of 
P and k 

-1 3/5 

P 



Pg = ump 



kBK{Mg 



(6) 



and assuming an ideal gas, 

ksT ^ k{Mgf'^P^I^, (7) 

keeping in mind K = K for the unmodified case. 
Furthermore , mot ivated in part by the simulations of 
iDolag et all (|2005D and lAscasibar et all (2003), we in- 
clude < 77 < 1 in ([5]) to allow for departures from strict 
hydrostatic equilibrium, in which the gas is supported 
by a combination of gas pressure and turbulent motions. 
The gr avitational potential of th e dark matter (DM) fol- 
lows a iNavarro. Frenk fc White (Il997| ) (NFW) profile, 
as motivated by N-Body si mulations, with a fixed con- 
centration parameter c = 5 (|Eke. Navarro, fc Steinmetd 



4 



l2001h . Finally, both the modified and unmodified models 
satisfy the same outer boundary condition, which is cho- 
sen to match the simulation results, but is related to the 
pressure required to resist the infall of baryonic matter 
at the viriol radius Pyir ~ hPDMV^ff, where /{, = ilb/rim 
is the cosmic baryon baryon fraction. 

The entropy floor is approximated by a shift in the 
entropy profile defined by: 



(8) 

As discussed in IVoit et all (|2002|) . this shift is a sim- 
ple phenomenological approximation of a uniform, high- 
redshift preheating. In particular it assumes that mass 
shells remain in the same order after the modification 
(this is reasonable due to that fact that entropy must 
be monotonically increasing for stability reasons), and 
that the effect of preheating is to increase the entropy by 
the fix ed, stated amount even after passing through the 
shock. IVoit et aTl (|2002[ ) examined this issue and found 
that this simple shift was a reasonable approximation (to 
within 20%) for smooth accretion. Of course, in a full 
cosmological simulation, the accretion is not generally 
smooth; however, as we will show, the shift approxima- 
tion works remarkably well. 

3.2. Matching the Unmodified Case 

For our comparison to be self-consistent, the most 
important consideration was our choice of unmodified 
distribution. This must match the results of the no- 
preheating simulation as closely as possible; for if our 
unmodified model does not match the simulations with- 
out preheating, we have no reason to expect that mod- 
ifications to that model will be consistent either. We 
therefore attempted the match the median entropy pro- 
file in two mass bins: low (4 < M200 < 6 x 10^^ h^Q 
Mq) and moderate (1 < M200 < 1-5 x lO" h:^^ Mq) 
mass clusters. All of the model predictions were made 
at the median mass in each bin; M200 = 5 x 10^'^ and 
-^200 = 1-3 X 10^'' h'^Q Mq in the low and moderate 
mass bins respectively. A summary of all the models 
considered is shown in Figure [H 

We first consider a set of ab-initio unmodified mod- 
els. Our fi r st cho ice was similar to the fidu cial model of 
IVoit et all (|2002l ) and lYounger et all (|2006 ^. in which the 
unmodified gas distribution is taken to follow the NFW 
dark matter density profile at the cosmic baryon density 
in strict hydrostatic equillibrium (77 = 1). Although it 
roughly reproduces the profile, in detail this NFW gas 
entropy profile was both too steep and incorrectly nor- 
malized, predicting a higher entropy for r ^ O.I-R200 than 
was seen in the simulations. It furthermore did not re- 
produce the core entropy seen in simulations, due to the 
divergent NFW density at small radii. Allowing r] to vary 
produced somewhat better agreement in the normaliza- 
tion and slope, but still predicted identically zero core 
entroy. 

We then tried a set of alternative models in an at- 
tempt to reproduce the core entropy see in the simu- 
lations : isothermal and a gas distribution with a flat 
core. Isoth ermal gas in hydrostatic equ ilibrium with 
b = 0.8 (see lMakino. Sasaki, fc Sutalll998f ) also had in- 
correct slope and normalization, with a core entropy that 
was too low. We then tried keeping the gas in hydrostatic 
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Fig. 2. — Unmodified entropy profiles from difl^erent models as 
compared to the median profile from the simulations (solid, thick) 
for clusters in both a low (left; 4 < A/200 < 6 X lO^^ h^^ Mq) 
and moderate (right; 1 < M200 < 1-5 X 10^* h^^ ^©) mass 
bin along with error bars representing the 1 — a variance about 
the median, and the associated fit from lVoit. Kay, fc BrvanI (120051) 
(solid, thin). The models shown; an NFW gas profile in hydro- 
static equilibrium (fg = 0.8, dotted), the isothermal model of 
IMakino. Sasaki, fc Sutol JToH, fg = 0.8 b = 0.8, dashed), gas with 
Pg ^ (1 + x)~^ {fg = 0.8, dot-dash), and an NFW gas profile in 
modified hydrostatic equilibrium (fg = 0.9 77 = 0.8, long-dash). 
The thick red line shows the fit to the simulation results used in 
our model. 

equilibrium, but modifying the gas distribution such that 
Pg ~ (1-1- x)~^. This too, did not fit the simulations, re- 
quiring a divergent temperature at the cluster center and 
a non-monotonically increasing entropy profile, both of 
which were unphysical. 

Since reproducing this median entropy profile from 
the no-preheating simulations with ab-initio models 
yielded no success, we instead chose to use a fit to 
the me dian entropy profil e in the AMR simulations of 
IVoit. K av. & Brvan'i l|2005[ )^. We parameterized the en- 
tropy profile in terms of fg, the enclosed gas fraction, 
with the following fit: 



Kifg)/K2m = 0.18 + 0.2/g + 1.5/2 



(9) 



where K200 = T2oo(/6200pcrit)"2/^/(Ai"^p) a nd T2Q0 = 
GM2o ofJ.m„/2r2oo, following the definitions in IVoit et al.l 
(|2002f ). We chose the outer boundary condition, once 
again based on the simulations, to be given by P200 = 

0- JT200.fbPDM(R20o)- 

IVoit. Kav. fc Brvaiil (|2005f ) showed that the entropy 
profiles in SPH and AMR simulations agree very well 
outside of 0.2r2oo while inside this radius, the SPH sim- 
ulations show steeper entropy profiles with lower (but 
non-zero) core entropies. The source of this discrepancy 
is not clear, al though the new entropy conserving version 
of Gadget2 (|Springe]|l2005f ) does show higher entropy 
cores (G. Yepes, private communication). In any case, 
the size of the difference at r = 0.1r2oo is quite small, 
about 25% and certainly not enough to resolve the dis- 
crepancy with observations discussed in section [H 

^ An alternative a pproach would have b een t o use polytropic fit s 
to simulations, as in lOstriker etal] | |2005D and IBode etaLl 1200^) . 
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Fig. 3. — Indicators of relative support from gas pressure (top) 
and turbulent motions (bottom) in the median profiles from the 
simulations for the same low (left) and moderate (right) mass bins 
from Figure [2] Shown are the results from four cases of injected 
entropy: no-preheating {Kq = 0; solid) and Kq = 78 (short-dash), 

155 (dot— dash) and 311 (long-dash) h^^^^ keV cm'^ . The error 
bars representing the variance about the median profile are, for 
clarity, only shown for the no-preheating case, but are similar for 
the other three. The thick solid line represent complete support 
against gravity, while the thick dashed line shows the relative con- 
tributions for J7 = 0.8. 

Finally, we note that our simulated clusters were not 
in strict hydrostatic equilibrium (see Fi gure [51), a result 
which is consistent wit h previous work (|Ascasibar et al.l 
l2003HDolag et al.ll2005l ). The resuh of this was that even 
with our fitted entropy profiles, we could not reproduce 
the simulated density and temperature profiles in detail 
for = 1.0. We found that setting 77 to 0.8, as suggested 
by the median profiles in Figure [3] (and previous work) 
generated a good fit to the gas profiles over a wide range 
of masses. 

3.3. Profiles 

Having found a reasonable entropy profile within our 
unmodified model, we then mo died these profile s ac- 
cording to th e pres cription in IVoit et all ()2002D and 
I Younger et al.l (|2006f) and computed the resulting model 
predictions for gas density, temperature, and entropy 
profiles in the same mass bins using matching entropy 

— 1/3 

shifts of = (no preheating), 78, 155, and 200 ft.yQ 
keV cm^ . In Figure|4l we compare the model predictions 
to the simulation results for the gas density pg, gas- 
weighted temperature Tg^, and gas-weighted entropy 
Kgw profiles. This comparison is made for the median 
profile in both the low and moderate mass bins. 

We find remarkable agreement between our analytic 
predictions and the simulation results over both mass 
bins for all three entropy shifts. The temperature profiles 
are underpredicted by the model, but never by more than 
~ 10% at i? ^ 0.1i?2oo> and never by more than ^ 20% 
out to -R200 ■ The entropy profile in the center is very well 
reproduced. There is also a tendency for the model to 
undcrprcdict the entropy at i? ^500 compared to the 
simulations, although again never by more than about 
20%. 
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Fig. 4. — Radial profiles at 2: = for the average gas density 
(pg), gas weighted temperature (Tg^), and gas weighted entropy 
(Kgnj) as a function of r / R200 as predicted by both the numeri- 
cal simulations (black, thick curve) and analytic model (red, thin 
curve) for four cases of injected entropy: no-preheating (Kq = 0; 
solid) and Kq = 78 (short-dash), 155 (dot— dash) and 311 (long- 
dash) h^Q^^ keV cm^ . The simulation results are the median 
profile for two representative mass ranges: 4 < M200 < 6 X 10^'' 
(left) and 1 < M200 < 1.5 X 10^* (right) h:^^ Mq. The analytic 
model predictions are taken at the median mass for each range in 
the numerical simulation: 5 X 10^'' (left) and 1.3 X 10^^ (right) h^Q 

Mq. 

We can therefore conclude, as one can see qualitatively 
in Figure [1] that preheating efficiently smoothes the ac- 
creted density distribution, making the accretion histo- 
ries of our simulated clusters more smooth. As a result, 
shifting the unmodified entropy profile is a good approx- 
imation to the effects of preheating in the simulations. 
Furthermore, we find that the model predictions and sim- 
ulation results converge as Kq increases. This occurs be- 
cause as Kg increases, the entropy input can effectively 
smooth larger and larger sub-haloes, and the simulations 
converge on the smooth accretion limit, making our mod- 
ified models a better and better approximation. 

3.4. The Intracluster Entropy 

From the comparison to the median profiles in Fig- 
ure [4l we find both qualitative and quantitative agree- 
ment between the analytic predictions and the simulation 
results. This is very encouraging, however comparing the 
median profiles is only a snapshot of the behavior of the 
intracluster entropy for the median cluster profile in two 
mass bins. The model predictions for the intracluster 
entropy over a much broader mass range contain a great 
deal of addition information. 

In Figure EJ we show the entropy-mass scalings for 
two characteristic entropy measures used in observations: 
the "core" entropy, taken at r = 0.1i?2oo (lower panel, 
ifo.ii?2oo)' and the "outer" entropy taken at r = i?5oo 
(upper panel, K^qq). For the core entropy, our analytic 
predictions are again in remarkable agreement with the 
simulations. However at high mass M200 ^ 3 x 10"'^^ h^Q 
Mq, our analytic predictions are somewhat lower than 
the simulation results. This behavior is slightly enhanced 
for higher values of Kq- Conversely, the models some- 
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M,„„ (h-„' 10" Mg) 

Fig. 5. — The intracluster gas entropy scalings at z = as 
predicted as predicted by both the numerical simulations and an- 
alytic model (points/ curves) . Shown are four different cases of in- 
jected entropy: no-preheating (Kq = 0; open triangles/solid) and 
K() = 78 (filled circles/short-dash), 155 (stars/dot-dash) and 311 

(filled triangles/long-dash) h^Q^^ keV cm^ . Show are the entropy 
at O.IR200 (Jfo.lfl2oo)' '"^'^ ^500 (-fsoo) as a function of M2oo- 




Fig. 6.— Same as Figure [5] for (-f^o.OlMg - -'^o.OiA/g )/-^^0 and 
(^O.SAfg - I<o.5Mg)/Ko, where A'o.OiAfg is the entropy at fixed 
gas fraction O.OlMg, -fCo.SMg is the entropy at fixed gas fraction 
O.bMg, K represents the no-preheating case, and Mg is the total 
gas mass at i?200. Note that a value of 1 represents no entropy 
amplification. 



what underpredict the outer entropy at all masses, with 
the exception of the very highest mass clusters where 
small numbers of objects make the trend difhcult to dis- 
cern. 

The agreement we see for the entropy of the preheated 
models is of particular interest and somewhat surprising. 
The effects of preheating on the entropy of simulated 
clusters in the cosmological simulation represents contri- 
butions from two competing effects. Entropy input into 
the ICM will tend to push out the accretion shock, which 
will in turn decrease the efficie ncy of he ating at the shock 
and decrease the core entropy (iBalog h. Babul, fc PattonI 
[T999HTozzi fc Norman|[200ll:lVoit et al 20031 However, 
as first suggested by IVoit et al.l ( 2003[ ). preheating will 
also tend to lower the density of accreted haloes. This 
may affect a partial transition from lumpy to smooth 
accretion that will actually increase the efficiency of en- 
tropy generation at the accretion shock, and with it the 
entropy of the system. Our simulations suggest that 
these competing effects tend to cancel each other out. 

In order to extract only the effects of entropy input 
on the simulations, we match the halos in each sim- 
ulation, and do a halo-by-halo comparison. For each 
cluster, we examine two ratios, or amplification factors: 

fcore = (-^^O.OlMg — Ko,oiMg)/Ko and f outer — (Ko.SMg — 

Ko.5m„)/Kq, where i^o-OiA/g is the gas-weighted average 



entropy (Tn 



-2/3 



at fixed gas fraction O.OIMg, Kq 



5M„ 



is the entropy at fixed gas fraction O.bMg, and Mg is 

the total enclosed gas mass at i?200- We use K to indi- 
cate the no-preheating case, and compute the ratio for 
each of the preheated simulations. These ratios represent 
the entropy offset between simulated clusters at fixed gas 
mass, in units of the entropy input Kq, at fixed gas mass 
in both the core and outskirts. The results of this com- 
parison are show in Figure [SI with fcore or f outer — 1 
representing no entropy amplification, and larger values 



of / implying more entropy amplification. 

We find only mild evidence of entropy amplification 
in Figure [51 The overall mean core and outer am- 

— 1/3 

plification factors for Kq — (78,155,311) hjQ keV 
cm^ are (/core) = (1.13,1.09,1.03) and {/outer) = 
(0.99, 1.03, 1.02) with associated variances acore = 
(0.36,0.36,0.18) and aouter = (0.34,0.21,0.16). The 
trend towards more amplification at lower masses for 
both ampl ification fact ors is qualitatively con sistent with 
IVoit et all (2003) and iBorgani et all (|2005l ). However, 
the mean amplitude of the amplification is < 15% for 
all values of Kq, and within a wide scatter is consistent 
with unity. Overall we find that on average, preheating 
does make entropy generation more efficient, but only by 
a small fraction. 

IBorgani et al.l ()2005f ) found amplification factors of 
~ 1.14 and 1.84 at R = 0.1 Rmr for a simulated group of 
mass 1.64 x 10^'^ h^Q Mq using entropy fioors of ~ 25 

and 100 h^Q^^ keV cm^ respectively. As noted above. 



we find smaller values for a much larger sample. One 
possible sou rce of this disc r epanc y is the nature of the 
simulations: 'Borga ni et al.l ()2005f ) include the effects of 
heating/cooling from a uniform UV background, star 
formation from a multiphase interstellar medium, and 
galactic winds powered by supernova feedback, whereas 
our simulations use a simple implementation of preheat- 
ing in which all of these processes are subsumed into 
the high redshift entropy modification. Another pos- 
sibility is the improved statistical power of our simu- 
lation s, which model more objects than IBorgani et al.l 
()2005D . Although there are groups of comparable mass 
in our sample with amplification fact ors that a re roughly 
consistent with the findings of Borg ajiTet al.l (2005), we 
also find that the mean amplification is within 1 — cr of 
unity, even for low-mass systems. However, perhaps the 
most likely explanation is the timing of the entropy in- 
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put.® In iBorgani et al.l ()2005D . the entropy was added 
much later, at z = 3, when many smaUer halos would 
have already formed. The gas in these halos is at high 
density and so, for a fixed amount of entropy, receives 
much more energy. This may have lead to large-scale 
outflows and shocking, which, in turn, may have caused 
the entropy amplification. In the simulations discussed 
he re, the input is done at rnuch ea rlier times {z = 10) as 
in iBialek. Evrard. fc MohJ (|2001h . because we want to 
isolate only the effects of preheating. 

We also find some evidence of turbulent mixing in the 
gas. This will will tend to equalize the entropy profile, 
raising the core and lowering the outer entropy. It has 
previously been noted that sim ulated clusters are no t 
in strict hydrostatic equilibrivim ( Ascasibar et al. 200^ , 
probably due to turbulence (Dolag et al. 2005). Figure [6] 
exhibits two telling trends: the core amplification fac- 
tor fcore tends towards unity with higher entropy input 
Kq, and the outer amplification factor f outer tends to- 
ward entropy suppression {fcore < 1) in systems with 
-^200 > 5 X 10^^ h^Q Mq. This is consistent with turbu- 
lent mixing, which will be less efficient at amplifying the 
core entropy at higher Kq as smoothing suppresses the 
turbulence, and more efficient at suppressing the outer 
entropy at higher mass as more of the substructure is re- 
solved and turbulence becomes more important. These 
two tends are, as we show in Figure [Sj borne out in our 
simulations; higher mass clusters with higher Kg tend to 
have more support against gravity contributed by turbu- 
lent motions in the gas. However, though we see trends 
that suggest turbulent mixing, wc concede that this is 
speculation and warrants a more in depth analysis, which 
we defer to future work. 

Therefore, our simulations show the combined effect of 
two processes that are difficult to disentangle. Classic 
entropy amplification will tend to increase both the core 
and outer entropy, and will tend towards higher fcore at 
lower mass. Turbulent mixing will suppress the outer en- 
tropy, tending towards lower f outer at higher mass, while 
simultaneously raising the core entropy more efficiently 
at higher Kq. We see all of these trends in our data. As 
a result, while the magnitude of the amplification favors 
the more subtle effects of turbulent mixing, we cannot 
firmly argue that either is dominant. Nevertheless, as 
Figure [6] makes clear, the size of either effect is not large. 

3.5. Scaling Relations 

In addition to the entropy scalings, we compare the 
model predictions to the simulation results for several 
observationally motivitated scaling relations. In Fig- 
ure [71 we present two scalings typically used by X- 
ray observers: the mass-temperature (for M200), and 
temperature-luminosity (for Lhoi, the bolometric X-ray 
luminosity) relations. In Figures [8] and [9] respectively, 
we present the mass and temperature scalings for two 
SZ observables defined above: the central decrement jjc 
and the FIADS SZ luminosity S^d\/f^. 

For the typical X-ray observables show in Figure [71 
we again find remarkable agreement between our model 
and the simulations. This is due largely to the agree- 
ment between the predicted and simulated core density, 
as the X-ray emission scales as ^ and is only weakly 

® We thank Mark Voit for pointing out this possibiUty. 
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Fig. 7. — Same as Figure O but for the temperature— luminosity 
and mass— temperature scahngs.The luminosity Z/f,oi is the bolo- 
metric X— ray luminosity and T is the emission weighted gas tem- 
perature. 

dependent on the temperature of the gas. Furthermore, 
this agreement suggests that the scalings in Figure [7| are 
largely insensitive to asymmetries in the simulations; we 
find that our assumption of spherical symmetry is a good 
approximation to the simulation results. 

In Figures [5] and [H we present the mass (left panel) 
and temperature (right panel) scaling of two typical SZ 
observables. Our analytic predictions are in good agree- 
ment with the simulation results for both scalings. We 
furthermore find that the central SZ decrement, in both 
the model and simulations, to be a far better indicator of 
cluster temperature than mass. We also confirm previous 
studies, which found the FIADS SZ integrated luminosity 
to be a robust, low-scatter mass and temperature indica- 
tor with roughly self-similar s caling (jda Silva et al.|[2000l : 
iMotl et al.l[2(mlNagaill200fih . This mass-FIADS SZ lu- 
minosity scaling is remarkably independent of the value 
of the entropy input Kq, while the temperature scaling 
exhibits weak Kq dependence. 

3.6. The Importance of the Core Entropy 

We have found very good agreement between our ana- 
lytic model predictions and the simulation results for the 
median density, temperature, and entropy profiles, the 
core and outer entropy scalings, and six different observ- 
able scalings. It is important to note, however, that this 
agreement came in large part by imposing an unmodified 
entropy profile derived from the simulations. None of our 
ab-initio unmodified models were in agreement with the 
simulation results. 

In particular, even when the slope and normalization 
of the entropy profile were matched closely, as in the 
case of an NFW gas profile with 77 — 0.8 (see Figure [2]), 
there was poor agreement with the simulation results. 
This was true even of the scaling of integrated quantities 
such as the Lboi — T and M200 — T relations, which were 
too steep and incorrectly normalized. At the same time, 
as the value of the entropy input Kq was increased, we 
found improved agreement with the simulations. 

This convergence is consistent with what we showed ex- 
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Fig. 8. — Same as Figure [71 . for the central SZ decrement j/c. 
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Fig. 9. — Same as Figure [T] for the frequency independent, an- 
gular diameter distance scaled integrated SZ luminosity Svd\/ fu. 

plicitly earlier; preheating efficiently smoothes the ICM, 
and therefore a simple shift of the unmodified entropy 
profile is a good approxi mation for the effects of pre- 
heating (jVoit et al.ll2003D . At the same time, an un- 
modified model with no core entropy was quantitatively 
and qualitatively at odds with simulation results, even 
when that unmodified entropy profile was a very good 
fit for r ^ 0.1i?2O0- And, though the value of this core 
entro py can vary depending on the numerica l technique 
used (|Voit. Kav. &: Brvan|[2005tlO'Shea etall feOOS). it is 
essential in setting the thermodynamic state and struc- 
tural properties of the ICM. The importance of this 
non-zero cor e entropy h as b een previously downplaye d 
(|Brvanll2000l : iVoit et al.|[200l iVoit. Kav. fc Brvanll200l , 
but here we find it to be a primary cause of disagreement 
between analytic model predictions and hydrodynami- 
cal simulations, even when the modification procedure is 
shown to be a good approximation. 



1 10 

T (keV) 

Fig. 10. — Temperature— luminosity and mass— temperature scal- 
ings at z = as predicted by both the numerical simulations with 
their associated power-law fits (black, points/curves). Shown are 
four different cases of injected entropy: no— preheating {Kq = 0; 
open triangles/solid) and Kq = 78 (filled circles/short-dash), 155 

— 1/3 

(stars/dot-dash) and 311 (filled triangles/long-dash) /i^q keV 
cm^ . Also plotted are the observations of IReiprich &: Bohringeil 
1I2OOI . red filled squares). 

4. COMPARISON TO OBSERVATIONS 

There has been some suggestion in the literature, mo- 
tivated by recent observations, that preheating mod- 
els which inv oke a univ ersal e ntropy floor are incor- 
rect in detail (|Ponman. Sanders on, k, Finoguenovl 120031: 
Pratt fc A rnaudI l2005t IPratt. Arnaud. fc P ointec outeaul 
2006). In particular, some authors have suggested that 
this brand of preheating predicts isentropic cores in low- 
mass systems that are not observed. In this section, we 
use the results of our simulations, which necessarily in- 
clude all the relevant dynamical and geometric effects^, 
as a test of the preheating scenario itself. Towards this 
end, we compare our results to observations of X-ray 
scalings (FigurefTU]) and both the core and outer entropy- 
temperature scalings (Figures [TTl [121 and [13]) . j 

In Figure 1101 we show the simulation results, along 
with best -fit power-laws, for the mass-temperature (up- 
per panel) and temperature-luminosity (lower panel) 
scaling relations as compar ed to the observations of 
IReiprich fc Bohringeil (|2002f ). As we mention above, 
the mass-temperature relation from the simulation is 
roughly independent of the value of the entropy in- 
put Kq. And, not surprisingly, with the exception of 
the no-preheating simulation all of the simulated mass- 
temperature relations are consistent with the observa- 
tional data. The temperature-luminosity relation, on the 
other hand, is very sensitive to the entropy input, as the 
bolometric luminosity scales roughly as LftoZ ^ Pg and 
preheating depresses and flattens the core density distri- 
bution. The temperature and luminosity are also both 
pure observables given a cosmological model - as opposed 
to M200 which Reiprich & Bohringcr ( 200^ estimate us- 
ing a n isothermal /?-model (jCavaliere fc Fusco-Femianol 
|1976[ ) - and therefore provide more robust constraints 

We note as a caveat that our simluations neglect radiative 
cooling, which may erode isentropic cores in low— mass clusters. 
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Fig. 11. — Same as Figure 1101 but for the outer clus- 
ter entropy K^qq. A lso show are the binned observations of 
IPonman. Sanderson. &: Finogucnov (2003, red filled squares). The 
solid line is the self-similar prediction from the simulations. The 
solid line is a power— law fit to the self-similar prediction from the 
simulations. 

on Kq. We find that an entropy input 155 ^ Kq ^ 311 

h^Q^^ keV cm^ wiU best fit the observational data. 

Figures [Til [HI and [13] show the outer (i^soo and 
^o.5i?.2oo) ^^"^ '^'^^'^ (^o.ii?,2oo) entropy-temperature scal- 
ing fro m the simul ations as compared to the observa- 
tions of iPonman. S anderson. & Finoguenov (2003) and 
iPratt. Arnaud. fc Pointecouteau (,2006i) . Our simulation 
results show the nearly isentropic cores in low-mass (or 
temperature) clusters: notice that the entropy is nearly 
mass-independent for low-mass clusters at the highest 
level of entropy input. In addition there is a break in 
the power-law scaling which shifts to higher mass with 
higher entropy input. 

We find that, while observations of the entropy at 
large radius in rich clusters are consistent our results, 
fully three-dimensional, non-radiative cosmological sim- 
ulations of the preheating scenario of cluster formation 
are inconsistent with observations in several ways. First, 
the entropy input scaling suggested by the temperature- 
luminosity relation is clearly inconsistent with observa- 
tions of the intracluster entropy; simulations of the best 
fit Kq range overproduce the core entropy and under- 
produce the outer entropy at fixed temperature. Second, 
the entropy-temperature scalings are themselves incon- 
sistent as compared to the observations; the best-fit en- 
tropy input value as determined by the core entropy un- 
derproduces the outer entropy, and vice-versa. Finally, 
we show the insentropic low-mass cores predicted by pre- 
heating as a function of Kq, which are not observed. 

As a caveat to our comparison to observations, 
we note that th ere has been s ome suggestion, 
most notably by iMcCarthv et all (|2004( l that ob- 
servational studies of the entropy profiles of clus- 
ters have been biased towards low core entropy sys- 
tems. This may, they argue, explain why some ob- 
servations appear to be inconsistent with high lev- 
els of e ntropy input (jPqnman. Sanderson, fc FinqguenoyI 
[200l IMcCarthv et all 120041: IPratt fc ArnaudI 120051 : 
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Fig. 12. — Same as Figure 1101 but for the outer cluster en- 
tropy -R'o.5i?2oo *° match the observations, which is similar to -fCsoo- 
Also s how are the observations of Pratt, Arnaud, & Pointccoutca^ 
II2OO6I . blue filled squares). The solid line is a power-law fit to the 
self-similar prediction from the simulations. The dashed line is the 
same r elation, approxim ately corrected t o the baryon fraction 
used in IPratt. Arnaud. fc Pointccoutea^ I I2006I) for comparison to 
the adiabatic scaling. We also note that, for three of the observed 
clusters, the measurement presented here is an extrapolation from 
~ 0.3-R200 and may be uncertain. 

IPratt. Arnaud. fc Pointecouteaul [20061 ) . However, while 
this is a possibility, we find dramatic and qualitative dis- 
agreement between the scalings in our simulations sim- 
ulations and observations which would be difficult to 
reproduce with a low-core entropy selection bias. We 
furthermore find disagreement at large radius (see Fig- 
ureslTTjandfT^. well outside the cool region. It is unlikely 
that a bias related to low entropy in the core would affect 
the outskirts. 

4.1. Including Radiative Cooling 

As an aside, we briefiy investigate the effects of ra- 
diative cooling. There has been some hope in the lit- 
erature that the shortcomings of the preheating sce- 
nario will be solved by a combination of preheating 
and radiative cooling (jVoit et al.l 120021 : IMcCarthv et al.l 
l2004f ). Unfortunately this is difficult to implement 
in a simulation because heating and cooling are prob- 
ably strongly tied togeth er (e.g., iKay et aT] 12004 

iMuanwong. Kav. fc Thomas! l2006t ). Nevertheless, we 
can imagine one simple and concrete extension to the 
straightforward preheating picture in which we include 
radiative cooling (but not star formation or feedback) at 
some late epoch, well after preheating has occurred. We 

therefore re-ran our Ko = 155 h^o^^ keV cm simula- 
tion ^ including the effe cts of radiative cooling for z < 3 
(see lO'Shea et al.ll2004D . We find that although the core 
entropy is modified, the entropy at -R500 is relatively in- 
sensitive to cooling and so the basic results above do not 
change. This is, however, another topic worthy of much 
more thorough investigation. 

5. CONCLUSION 

We present high-resolution hydrodynamic simulations 
of the preheating scenario of cluster formation, in which 
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Fig. 13. — Same as Figure 1111 but for the core clus- 
ter entropy ^^0.1^,00- Also show arc the binned observations 
of Ponrnan. Sanderson, & FinoKuenov (2003, red filled squares) 
and Pr att. Arnaud. fc Point ccouteau (2006t, blue filled squares). 
The solid line is a power-law fit to the self— similar predic- 
tion from the simulations. The dashed line is the same re- 
la tion, approximately corrected to the b aryon fraction /(, used 
in IPratt. Arnaud. &: Pointecouteaul II2006I) for comparison to the 
adiabatic scaling. The offset between the two o bservations 
is like ly due to the fact that IPratt. Arnaud. &: Point ecouteaul 
II2OO6I) resolve the temperature profiles of their clusters, 
whereas IPonman. Sanderson, fc Finoguenovl (120031) do not. This 
may lead to an over-estimatio n of the core entropy by 
IPonman, Sanderson, &: Finoguenovl (120031 ) 

the effect of SN and/or AGN heating is approximated as 
a universal increase of the entropy, implemented at high 
redshift. We focus on two aspects: (1) how well analytic 
models of entropy input reproduce the simulations re- 
sults, and (2) how well the preheating simulations match 
current cluster observations. 

We fi nd that analytic mode ls of pr eheating following 
those of IVoit et all (20021 ) and lounger et al. (200^ are 
remarkably successful at reproducing the results of hy- 
drodynamical simulations in detail. This agreement ex- 
tends from scalings of integrated quantities such Lhoi — T, 
M200 — T, and SZ scalings for j/c and the FIADS in- 
tegrated SZ luminosity, to core and outer entropy scal- 
ings, and finally to the median entropy, temperature, and 
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density profiles of both low and moderate mass clusters. 
However, this agreement is dependent on the right choice 
of an unmodified profile; in particular one with non-zero 
core entropy. 

When wc look in detail at the simulation results, we 
do not find strong evidence for entropy amplification, the 
process in whi ch the entropy is amplified even beyond the 
level inj ected (IVoit e t al."2G03f). This differs from the re- 
sults of iBorgani et al. (2005), who looked at a small set 
of clusters and groups. The reasons for this are not en- 
tirely clear although their preheating model differed from 
that used here (in particular the entropy input occured 
at a lower redshift). 

Our simulations also show that the preheat- 
ing scenario for cluster formation is in confiict 
with recen t observations of the X-ray scalings 
(jReiprich fc Bohringer 200 1) and intraclu ster en- 
tropy ~(iPonnian. Sanderson, fc Finoguenovl 120031: 
Pratt fc ArnaudI 120051 : IPratt. Arnaud. &: Pointecouteaul 
20061 ). The best-fit value of the entropy input Ko, 
determined from the X-ray scalings (see Figure ITU)) , 
underproduces the entropy at i?5oo and overproduces the 
entropy at 0.1i?2oo at fixed temperature (see Figures [TT] 
and [T5)) . No single value for the entropy input due 
to preheating can simultaneously match the observed 
entropy values at both the core and the outskirts. This 
result is in qualitat ive agreement with previous work 
(jBorgani et al.l [20051 ) . The inclusion of radiative cooling 
at late times (without any additional feedback) does not 
appear to change this result. Therefore, we find this 
simple preheating scenario to be in disagreement with 
observations, implying that other, more sophisticated 
treatments are needed. 
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